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Abstract 

We explore the natural limit of binomial reducibility in nuclear multifragmen- 
tation by constructing excitation functions for intermediate mass fragments 
(IMF) of a given element Z. The resulting multiplicity distributions for each 
window of transverse energy are Poissonian. Thermal scaling is observed in 
the linear Arrhenius plots made from the average multiplicity of each element. 
"Emission barriers" are extracted from the slopes of the Arrhenius plots and 
their possible origin is discussed. 
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Emission of multiple intermediate mass fragments (IMF), 3 < Z < 20, is an important 
decay mode in heavy-ion collisions between 20A and lOOA MeV Despite extensive 

studies, the nature of the fragmentation process, whether statistical or dynamical, remains 
an open problem. An historic overview of low energy reactions shows that the emission prob- 
abilities and excitation functions are by far the best observables in distinguishing between 
statistical processes (dominated by phase space, as in the case of light particle evaporation 
and fission), and prompt, dynamical processes (like direct reactions) 0. Indeed, several 
aspects of nuclear multifragmentation may be understood in terms of nearly independent 
fragment emission from multifragmenting sources with thermal-like probabilities f^-^]. 

It was found 0-0] that the experimental Z integrated fragment multiplicity distributions 

are binomially distributed in each transverse energy (Et) window, where n is the number 
of emitted fragments and m is the number of throws. The transverse energy Ef is calculated 
from the kinetic energies Ei of all the charged particles in an event and their polar angles 
9i, as Et = J2i Eisiia^ 9i. The extracted one-fragment emission probabilities p give linear 
Arrhenius plots (i.e. excitation functions) when log 1/p is plotted vs 1/y/E^. If the excitation 
energy E* is proportional to Et and consequently, the temperature T to ^/E^, these linear 
Arrhenius plots suggest that p has the Boltzmann form p oc exp(— 5/T) 

Similarly, the charge distributions for each fragment multiplicity n are observed to be 
reducible to a single charge distribution and to be thermally scalable Also, the experi- 
mental particle-particle angular correlation is reducible to the individual fragment statistical 
angular distributions and thermally scalable 0. 

The appeal of this comprehensive picture is marred by a number of open problems. 
One problem, which will be dealt with here, is that the binomial decomposition has been 
performed on the Z-integrated fragment multiplicities, typically associated with 3 < Z < 20. 
Thus, the Arrhenius plot generated with the resulting one fragment probability p is an 
average over a range of Z values. Fortunately, it has been shown that the Arrhenius plots 
should survive such a Z averaging, and yield an effective "barrier" (slope) dominated by 
the lowest Z value 0-^]. However, this procedure clearly implies a substantial loss of 
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information, and renders the binomial parameters p and m difficult to interpret. 

In light of the above considerations, an analysis of the multiplicities for each fragment 
Z value may solve many of these difficulties. Furthermore, it has been pointed out that a 
binomial distribution could be distorted by the averaging associated with the transformation 
E* Et leading to possibly incorrect values of m and p [|TU|. However, it has been shown 



that, while p and m separately can conceivably be distorted by the transformation, the 
average multiplicity (n) = mp is far more resistant to the averaging process It 
would be useful if a way could be found of avoiding the individual extraction of p and m 
while retaining the possibility of constructing an Arrhenius plot. 

In this letter, we analyze the experimental fragment multiplicity distributions for each 
individual fragment Z value. We show that they are Poissonian. The associated mean mul- 
tiplicities for each Z give linear Arrhenius plots from which the corresponding Z dependent 
barriers can be extracted. The physical dependence of these barriers on Z may shed light on 
the fundamental physics associated with multifragmentation, as fission barriers have done 
for the fission process. 

The effect of restricting the fragment definition to a single Z value is rather dramatic. 
In Fig. ratios of the variance to the mean as a function of Et are given for a number of 
Z values, and for the case Z > 3 P]. For individual Z values the ratios are very close to 
one, while for the Z integrated case there is a sagging at large Et. The explanation for these 
features can be found by recalling that for a binomial distribution 

(n) = mp; = mp{l — p); —— = 1 — p. (1) 

{n) 

For p — > 0, the ratio cr^/ (n) 1. This is the Poisson limit. When extensive summation 
over Z is carried out, the elementary probability p increases sufficiently at the highest values 
of Et so that the Poisson distribution is replaced by the more general binomial distribution. 
On the other hand, the restriction to any given Z value decreases the elementary probability 
p so dramatically that the above ratio effectively remains unity at all values of Et and the 
distributions become Poissonian: 



where (n) is (ra) (-E*). We show the quahty of the Poisson fits to the multiphcity distribution 
in Fig. ^ These Poisson fits are excellent for all Z values starting from Z=3 up to Z=14 over 
the entire range of Ef and for all the reactions which we have studied. Thus we conclude 
that reducibility (we should call it now Poissonian reducibility) is verified at the level of 
individual Z values for many different systems. Incidentally, for the Xe induced reactions, 
there is an excellent overlap of the data sets for different targets as a function of Et. They 
all follow the Poisson fit to the Au target data. The probabilities Pn and the range of Et 
increase with the increasing target mass from V to Au, as they must if Et is a reasonable 
measure of the dissipated energy. 

The experimental observation of Poissonian reducibility directly implies that IMF pro- 
duction is dominated by a stochastic process. Of course stochasticity falls directly in the 
realm of statistical decay, either sequential or simultaneous (see section 5.5 of ref. or 
ref. 0]). It is less clear how it would fare within the framework of a dynamical model. 

In order to verify thermal scaling i.e. if the emission probabilities are thermal, we generate 
Arrhenius plots by plotting log (n) vs 1/ y/E^. Here, as in previous works, we assume that 
Et oc E* and that E* oc T^, according to the simplest strongly degenerate Fermi gas 
dependence at constant volume. We are of course aware that high excitation energies and/or 
lower densities can lead to deviations, which may well be looked for in the future. We 
expect (n), like p, to be of the form (n) = F{T, ...)e~^^^, where the specific form of the pre- 
exponential factor depends ultimately on whether a reaction theory or a chemical equilibrium 
description will prevail. We use the Arrhenius plot in the traditional spirit of evidentiating 
the leading T dependence contained in the exponential. The top four panels of Fig. ^ gives 
a family of these plots for four different reactions. Each family contains Z values extending 
from Z=3 to Z=1A. The observed Arrhenius plots are strikingly linear, and their slopes 
increase smoothly with increasing Z value. One slight exception is the large Z (> 10) data 
for Xe+Cu. At high Et, the data deviates from the linear dependence observed elsewhere. 
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For this smaller system, it is conceivable that charge conservation constraints lead to this 
behaviour. The overall linear trend demonstrates that thermal scaling is also present when 
individual fragments of a specific Z are considered. Even apart from the linearity of the 
Arrhenius plots, important information is already contained in the range covered by the yield 
of individual fragments over the range of Ei shown in Fig. ^. For processes not dominated by 
phase space (e.g. low energy direct reactions), one expects the excitation function to depend 
weakly upon excitation energy. Typically the cross sections vary by factors of a few. In the 
present data, the mean multiplicity (n) varies with Et by one to two orders of magnitudes. 
This is a strong evidence for the involvement of the internal degrees of freedom typical of 
high barrier statistical decays. 

The advantage of considering individual Z selected fragments is readily apparent. For 
any given reaction, both Poissonian reducibility and thermal scaling are verifiable not just 
once, as in the binomial analysis, but for as many atomic numbers as are experimentally 
accessible. Take for example the Ar+Au reaction (i?/y4=110MeV) shown in the top right 
panel of Fig. ^. For this specific reaction, we can verify both reducibility and thermal 
scaling for 12 individual atomic numbers. Since there are 29 Ei bins, Poissonian reducibility 
is tested 29 times for each Z value, i.e., 12 x 29 = 348 times for this reaction alone. Including 
all the cases shown in Fig. |^, we have tested Poissonian reducibility 936 times. This is an 
extraordinary level of verification of the empirical reducibility and thermal scaling with the 
variable Et- 

Two added bonuses arise from this procedure. 

1) The criticism has been raised that the linearity of the Arrhenius plots arises from 
an autocorrelation, since the complex fragments also contribute to Et [|1^]. In the present 
analysis this criticism can be dismissed, since each individual Z contributes a vanishingly 
small amount to Et (< 5%), even in the region of maximum yields. Still, to be sure that 
there is no autocorrelation in Fig. |^, we have repeated the analysis, for Xe+Au at 50A MeV, 
by: i) removing from Et the contribution of the individual Z (E^) that we have selected 
(Fig. H, bottom left panel); ii) by using only the Et of the light charge particles, Ef"'-^^ 
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(Fig. 0, bottom right panel) . In both cases, the Arrhenius plots remain linear over almost 
the entire range of Ef and cover 1 to 2 orders of magnitude. Quantitatively, the rate of 
change of the slopes with Z remains the same regardless of the definition of Et, as shown in 
the top panel of Fig. ^. This behaviour is expected if the slopes are related to some physical 
barriers. 

In our attempt to avoid autocorrelation by excluding from Et all IMFs {E^^^) or the Z 
value under investigation {E^), we have introduced another kind of distortion. Excluding 
from Et all fragments of charge Z to produce Et necessarily requires that for those events 
where E^ ^ Et, the yield nz — > 0. This produces the visible turnover of the Arrhenius 
plots in the bottom panels of Fig. ^ (the same argument also applies to E^'^^). It has been 
verified experimentally that the maximum values of the new Et scale do indeed correspond 
to events in which the contribution from a given Z (or all IMFs) is absent. 

2) The extracted elementary probability is now (n) = (mp) which, contrary to p and m. 



is very resilient to any averaging associated with the transformation from E* to Et [p!0| , p |. 

It may be worth reminding the reader that this procedure does not contradict binomial 
reducibility. To the contrary, it represents its natural limit for small values of p, and it 
expands its applicability by considering each Z value individually. In going from binomial 
to Poissonian distributions, the price one pays is the loss of the parameter m. While in many 
ways this is a convenient result, it actually implies a loss of scale. In the time sequential 
interpretation of multifragmentation |^ this implies a loss of information about the time 
window during which multifragmentation occurs in units of the natural channel period, or 
the unit time to which the elementary probability is referred. In the spacial interpretation, 
one loses information about the total mass of the source 0. 

Poissonian reducibility and thermal scaling do not contradict recent observations regard- 
ing the role of reaction dynamics in the formation of the hot primary sources p^-p!9|. In 



particular, the experimental scaling is not affected by the presence of multiple sources 
and the analysis presented here is a powerful test to establish the degree of thermalization 
in the late stage of the reaction. Kinematic variables seem to retain spatial-temporal in- 
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formation about the reaction dynamics |lT5|-pO| while the associated emission probabihties 



seem to demonstrate, as verified nearly a 1000 times in the present work, the role of phase 
space in describing the decay of the sources. 

Returning to the Arrhenius plots for individual atomic numbers, it is straightforward 
to obtain the values of the slopes from Fig. ^ as a function of Z . The interpretation of 
these slopes as "emission barriers" is very tempting. If we had the correct excitation energy, 
rather than Et^ we could obtain the actual barriers as a function of Z. Unfortunately we 
are limited to our running variable E't, and to the assumption of its proportionality to E* . 
However, the many linear Arrhenius plots shown here cannot be easily explained without 
invoking this proportionality. Therefore, with the necessary caution, we explore the possible 
meaning of these "barriers". A plot of these barriers as a function of Z is potentially rich 
in information. The extracted barriers are shown in Fig. ^ (bottom panel). The barriers 
appear to increase linearly with Z at low Z values and tend to sag below linearity at higher 
ones. 

One could wonder about the role of surface energy on the origin of these barriers. Frag- 
ments might be thought as forming by coalescence into a relatively cold and dense nuclear 
drop out of a hot diluted source. The appearance of a substantial surface energy for the frag- 
ment would suggest barriers proportional to Z^/^ (A^/^). If this were true, then one would 
expect the barrier for each Z to be nearly independent of the system studied. Unfortunately, 
since the relation between Et and excitation energy is unknown, the absolute values of our 
barriers are also unknown. By normalizing all systems at Z=6 and using the Xe+Au at 
E/A=50 MeV as the reference, one observes barriers that are indeed fairly independent of 
the system (Fig. H, bottom panel). Another possibihty is to compare the dependence of the 
barriers on Z to that of the conditional barriers measured at low energy |^ (black dots). 
Their similarity with the multifragmentation barriers is dramatically illustrated. While the 
Coulomb-like Z dependence of these barriers is suggestive, we should remark that these are 
emission barriers rather than Coulomb barriers. Thus the dominance of the Coulomb term 
is by no means obvious. 



In conclusion, Poissonian reducibility and thermal scaling of individual fragments of a 
given Z have been observed experimentally for several different systems at bombarding 
energies ranging from 50 to 110 MeV/nucleon. The high level of verification strongly sup- 
ports the stochastic/statistical nature of fragment production and provides a clear signal for 
source(s) thermalization in the late stage of the reaction. Slope parameters were extracted 
from the Arrhenius plots. The interpretation of these slopes as emission barriers, originating 
either from Coulomb or surface terms, or both, still needs to be explored. If the physical 
significance of these Z dependent "barriers" must remain lamentably open, there is at least 
the distinct possibility that important physical information is contained therein. Data with 
isotopically resolved hght charged particles and IMFs are needed to further investigate these 
phenomena. 
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FIGURES 



FIG. 1. The ratio of the variance to the mean number of Li, C, O and Ne fragments (sohd and 
open symbols) emitted from the reaction '^^Ar+^^'^Au at E/A=1W MeV. The star symbols show 
the same ratio for ah IMFs (3 < Z < 20). 



FIG. 2. The excitation functions P„ for carbon (left column) and neon emission (right column) 
from the reactions ^SAr+^^^Au at £^7^=110 MeV (top panels) and ^'^^Xe+^'^V,''''^Cu,^^Y,^^'^ Au 
(bottom panels). The lines are Poisson fits to the gold target data. 



FIG. 3. Middle and upper panels: The average yield per event of different elements (symbols) 
as a function of Bottom panels: The Xe+Au data at 50A MeV are replotted using the 

transverse energy of all charged particles excluding the Z that we have selected, (left), and 
(right) that only of the light charged particles, E^^^ . The lines are fits to the data using a 
Boltzmann form for {nz)- 



FIG. 4. Top panel: Slopes of the Arrhenius plots, normalized to Z = 6, for Xe+Au at 50A 
MeV as a function of Z using the indicated definitions of Et. Bottom panel: The Z dependent 
"barriers" (the slopes of the Arrhenius plots in Fig. |3|). The "barriers" have been scaled relative to 



Z = 6 of the Xe+Au data. Black dots are low energy conditional barriers from ref. |21| normalized 
to Z = 6 of Xe+Au. 
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